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ABSTRACT 

We measure the geometry, topology and morphology of the superclusters in mock 
SDSS catalogues prepared and reported by Cole et al.(1998). The mock catalogues 
refer to rCDM and ACDM flat cosmological models and are populated by galaxies 
so that these act as biased tracers of mass, conforming with the observed two-point 
correlation function measured using APM catalogue on scales between 1 to 10 h~^Mpc. 
We compute the Minkowski Functionals (hereafter, MFs) for the cosmic density fields 
using SURFGEN (Sheth et al.2003) and use the available 10 realizations of tCDM to 
study the effect of cosmic variance in estimation of MFs and Shapefinders; the statistics 
derived from MFs, and used to study the sizes and shapes of the superclusters. The 
MFs and Shapefinders are found to be extremely well constrained statistics, useful in 
assessing the effect of higher order correlation functions on the clustering of galaxy- 
distribution. We show that though all the mock catalogues of galaxies have the same 
two-point correlation function and similar clustering amplitude, the global MFs due 
to rCDM show systematically lower amplitude compared to those due to ACDM, an 
indirect, but detectable effect due to nonzero, higher order correlation functions. This 
enables us to successfully distinguish the two models of structure formation. 

We further measure the characteristic thickness (T), breadth (B) and length (L) 
of the superclusters using the available 10 realizations of rCDM. While T^B and T, 
Be [1,17] h~^Mpc, we find the top 10 superclusters to be as long as 90 h~^Mpc, with 
the longest superclusters identified at percolation to be rare objects with their length 
as large as 150 h~^Mpc. The dominant morphology of the large superclusters is found 
to be filamentary. The thickness, breadth and planarity of the superclusters follow well- 
defined distributions which are different for the two models. Thus, these are found to 
be sensitive to the cosmological parameter-set and are noted to be candidate statis- 
tics which can compare the rival models of structure formation. Further, the longest 
structures of rCDM are found to be significantly longer than those in ACDM. Finally, 
mass and volume-weighted dimensionless Shapefinders - Planarity and Filamentarity 
- are found to be well-constrained statistics useful to discriminate the two models. 
We note some interesting effects of bias and stress the importance of incorporating 
realistic treatment of bias in preparing and analysing the mock catalogues. 
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1 INTRODUCTION 

The present era marks a remarkable milestone in the history 
of Cosmology due to an alround increase in the extragalac- 
tic database by several orders of magnitude over last few 
years. These data have improved both in quantity and in 
quality, due to which theoretical models of the Universe are 
directly confronted with the observations, either to be vali- 
dated or to get more refined. The primordial density fluctu- 
ations as revealed imprinted on the otherwise smooth and 



isotropic microwave background radiation (MBR) have been 
measure d with unprecedente d accuracy and resolution by 
WMAP iBennett et al.l2003h . thus enabhng us to obtain the 
most reliable estimates of the cosmological parameters and 
confirming that the geometry of the Universe is indeed flat. 
An independent constraint from the study of extragalactic 
distances using type-la supernovae seems to indicate that 
we live in a low matter-density, cosmological constant dom- 
inated Universe, with Qa «0.7. A variety of observations 
like, the power spectrum study of the MBR, the clustering 
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properties of the galaxy-distribution and the rotational ve- 
locity measurement s in normal-size gal axies as well as in the 
clusters of galaxies (|Kneib et alJl200^ . all seem to indicate 
that most of the remaining fraction of the matter-density 
is in the form of dark matter which is nonrelativistic, i.e., 
cold in nature. Thus, the cosmological constant dominated. 
Cold Dark Matter model (ACDM) seems to be the best 
fitting cosmological model of the Universe that we inhabit. 
However, several anamolies have surfaced which cannot be 
easily reconciled with the hypothesis of cold dark matter. 
For example, an order-of-magnitude lower than predicted 
number density of dwarf or low surface brightness galaxies 
observed in the local Universe, the mismatch between the 
predicted and the observed density profiles of these galax- 
ies, the lack of sufficient substructure in normal size galaxies, 
going against the standard hierarchical model of structure 
formation, etc. are not readily explained within the standard 
CDM paradigm and secondary physical processes are sought 
to explain many of these anomalies. Hence, the issues related 
to structure formation and background cosmology are likely 
to have a vigorous interplay with the physics of galaxy for- 
mation, a process which may lead to deeper insights into 
our ideas about the cosmological model(s) of the Universe. 
It suffices, hence, to devise algorithms and statistics which 
enable us to compare our theoretical predictions with a va- 
riety of observations of the Universe. 

The developments on studying the microwave back- 
ground radiation are further complemented by properties of 
the large scale structure revealed by increasingly deep and 
voluminous redshift surveys. Indeed, redshift surveys con- 
tinue to occupy a central place in our efforts to understand 
the Universe for, as opposed to CMBR which gives us a view 
of the early Universe, these provide us with a view of the 
matter-distribution in our immediate extragalactic neigh- 
bourhood. In addition to revealing the large-scale distribu- 
tion of galaxies, a redshift survey provides us with a statis- 
tically fair sample to estimate many properties associated 
with various galaxy-populations, e.g., the lumin osity func- 
tion llBlanton et. aljbooal : fNorberg et al."2002aV the clus- 
tering of galaxies of various "types" |NorbcrE ct al. 2002lJ, 
the mass-to- light (M/L) ratio in galaxies of different mor- 
phological or spectral types (for example, Brainerd & Spe- 
cian 2003), etc. Further, with the help of efficient algorithms, 
it is possible to identify groups and clusters of galaxies and 
assess the role which their environments play in the evolu- 
tion of their constituent galaxies ^. Such findings provide 
crucial feedback to our efforts in understanding the physics 
of galaxy formation. For a detailed review highlighting po- 
tential applications and importance of redshift surveys in 
the study of the Un iverse, see a recent review-article by 
iLahav fc Su'tolj2003D . 

Whereas the paradigm of structure formation due to 
gravitational instability has been explored to a great extent, 
both analytically and using N-body simulations, these stud- 
ies have mainly focussed on the dynamics of the dark matter. 
However, we now have richer quality and quantity of data 
with which to confront our theoretical models. For example, 
the 2 degree field galaxy redshift survey (2dFGRS) contain- 

^ See lEinasto et al.l l2003l) for a study of the clusters and super- 
clusters in the SDSS Early Data Release. 



ing about 250 000 galaxies and covering about 1500 square 
degree of the sky is complete (Colless et al.2003), while the 
Sloan Digital Sky Survey (SDSS) team has recently made 
its first data relea se (FDR) reporting the redshifts of about 
180 000 galaxies llAbazaiian et al.ll2003ll . When complete, 
the SDSS will cover about a quarter of the full sky and re- 
veal the distribution of a million galaxies upto 600 h~^ Mpc, 
a depth comparable to that attained by 2dFGRS. 

In light of the above developments, the attention of 
the cosmology community has steadily shifted to simulating 
the formation of galaxies in dynamically evolving dark mat- 
ter background. In this context, it is worth mentioning two 
widely adopted a pproaches - the semi-analytic galaxy for- 
mation methods lcol^^L|j200^jBensoii et al J 120021) and 
SPH methods llWeinberg. Hernauist fc Katzl bOoJT How- 
ever, these two approaches do have their own limitations 
(Berlind et al.2003), and important insights have been ob- 
tained by a complementary approach based on the Condi- 
tional Luminosity Function (see Yang, Mo, van den Bosch 
& Jing 2003 and a series of earlier papers by this group). 

The catalogues of "galaxies" prepared using the meth- 
ods advocated by these teams are to be tested against the 
large body of observations due to deep, upcoming redshift 
surveys. Such a comparison should have two main aims: 

(i) To employ various population statistics which reflect 
the structural properties of the ensemble of galaxies. These 
are mainly the kinematic measures defined for a given sam- 
ple of a large number of galaxies. The luminosity function, 
the mass-to-light ratios of galaxies of various morphological 
or spectral types, mass-to- light ratios of groups/clusters, the 
morphology-density relation etc. form this class of measures 
and are chiefly used to validate/invalidate a mock catalogue. 

(ii) To assess the large-scale distribution of galaxies and 
to quantify the clustering of the ensemble of galaxies. The 
two-point correlation function, the three-point correlation 
function, the pair-wise velocity dispersion, etc. are the most 
widely used statistics belonging to this category. These 
statistics reflect the dynamical history and clustering prop- 
erties of the matter-distribution. 

The list of statistics which fall in the latter category is, 
however, not complete. The statistical properties of the pri- 
mordial gaussian random phase density fleld are fully quan- 
tifled in terms of its power spectrum, or its Fourier counter- 
part, the two-point correlation function. However, a subse- 
quent gravitational dynamics introduces phase correlations 
among the neighbouring density modes of this density field 
(Sahni & Coles, 1995). This is accompanied by a transfer of 
powe r from larger scales to smaller scales (iHamilton et alJ 
I1991I) as a result of which the entire hierarchy of correlation 
functions becomes nonzero. Hence, a full quantification of 
the large scale cosmic density field at today's epoch would 
ideally require us to have knowledge of all the higher or- 
der correlation functions. Under appropriate assumptions, 
analytic treatment goes as far as predicting the effect of 
gravitational dynamics on the 3-point correlation function, 
or at best the 4-point correlation function, but not beyond. 
It is to be noted that even this treatment is within the per- 
turbative regime and breaks down when the density con- 
trast exceeds unity. However, we note that a mock cata- 
logue could not be satisfactorily compared with the real 
cosmic density field merely by a subset of correlation func- 
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tions for which analytic predictions are available. In fact, 
most of the N-body simulations are structure-normalised 
so as to reproduce the current rich cluster abundance ob- 
served in our neighbourhood. In addition, mock catalogues 
constructed from N— body simulations are usually required 
to reproduce the observed two-point galaxy-Kalaxy cor- 
relation function dCole. Hatton. Weinberg, fc FrenMllQQSt 
lYang. Mo. van den Bosch fc Jind |2003^ . Thus, the mock 
catalogues of various cosmological model(s) agree with the 
large scale structure (hereafter, LSS) of the Universe by con- 
struction, and cannot be discriminated from one another 
or from the observed LSS of the Universe on the sole ba- 
sis of the 2-point correlation function. We might add that 
higher order correlation functions (beyond, say, the three- 
point correlation function) are cumbersome to estimate for 
a system comprising of large number of particles. Further- 
more, the correlation functions offer us little intuitive insight 
to guide us in interpreting the results. Thus, distinctly differ- 
ent statistics are needed to give us an handle on the effect of 
higher order correlation functions. It is desirable that these 
statistics simultaneously offer us an intuitively clear inter- 
pretation as well as allow us for their numerically robust 
estimation. 

In this paper, we propose to use one such class 
of me asures - the Minkowski Functionals (hereafter , 
MFs) jMecke. Buchert fc WagneJ Il994l: ISheth et alJl2003l) 
and related morphological statistics, the Shapefinders 
iSahni. Sathvaprakash fc Shandarinll99^ - to quantify and 
characterise the large scale distribution of galaxies^. Our 
main motivation in using the MFs for this purpose is that 
(1) MFs have been shown to be dep endent on the entire hi- 
erarchy of correlation functions (Me cke. Buchert fc Wagneil 
ll994l:ISchmalzing|l999ll . and (2) the MFs in both 2-D and in 
3— D have physically important interpretation, namely that 
out of (n+l) MFs available in n— D, the first n of them char- 
acterise the geometry of the excursion set under question 
(these could be superclusters or voids of the cosmic web, for 
example) and the remaining MF is sensitive to the topology 
or the connectedness of the associated hypersurface. In re- 
cent times, there is even more motivation for employing MFs 
for the quantification of LSS after Matsubara (2003) derived 
the MFs for cosmic density fields in the weakly nonlinear 
regime using 2°'* order Euler perturbative expansion. With 
the availability of large datasets like 2dFGRS and SDSS, 
it should be possible to test these predictions by working 
with the cosmic density fields smoothed at sufficiently large 
length-scales ^. 

Various attempts have been made to quantify the cos- 
mic web using MFs iBharadwai. et. aL I I2OOOI : ISheth et alJ 
[200 3; ?). The efficacy of MFs in quantifying cosmic den- 
sity fields as well as in making a morphological survey of 
the structural elements of the cosmic web - superclusters 

^ See also, Doroshkevich et al.2003 for a morphological study of 
the SDSS (EDR) using minimal spanning trees. 
^ In fact, one could attack long lasting questions related to bias 
on large length-scales (of about a few tens of Mpc) employing a 
completely different line of approach of estimating MFs for cosmic 
density fields smoothed on, say a few ten's of Mpc and comparing 
the results with their theoretical predictions for the dark matter. 
Sec Hika gc. Taruva fc Suto ( 2003) for an application involving 
genus. 



- has thereby been abundantly demonstrated. Recently, 
iHikaee et alJ (jioO^ p rovided M i nkow ski Functionals for 
SDSS (EDR). Further. iBasilakosI ll2003t) used Shapefinders 
to study the morphology of the superclusters in SDSS(EDR) 
and found the dominant morphology of the superclusters 
to be filamentary. In the mean time, important advances 
have been made concerning techniques to calculate the MFs. 
The grid-based techniques ba sed on the Crofton's formulae 
or the Koenderink invariants JSchmalzing fc Buchertlll99iD 
have been recently complemented by a powerful, triangula- 
tion based surface-generating method which self-consistently 
generates isodensity contours for a density field defined on 
a grid. This method also makes possible an online computa- 
tion of the MFs. The method, the associated software SUR- 
FGEN, its accuracy and application to a class of N— body 
simula tions were reported in a recent paper bv lSheth et alJ 
J2003h . 

In this paper, we analyse the LSS in the SDSS mock cat- 
alogues corresponding to two models of structure formation 

- ACDM and rCDM. These mock catalogues were prepared 
by Cole et al.(1998) (hereafter, CHWF) based on various 
physically motivated biasing schemes represented in terms 
of analytic functions. The biasing schemes were used to pre- 
scribe the "formation-sites" of the galaxies for a given cos- 
mological density field. Despite being superseded by semi- 
analytic and Conditional Luminosity Function based tech- 
niques to develop more realistic mock catalogues, till date 
these catalogues have several advantages to their credit. This 
is because, these catalogues fulfill several kinematic as well 
as dynamical constraints. For example, the survey geometry 
and the hj band luminosity functions for both 2dFGRS and 
SDSS have been incorporated, ensuring an accurate repro- 
duction of the expected radial selection function in both the 
cases *. In addition, these catalogues mimic the survey geom- 
etry of both 2dFGRS and SDSS, thus pausing the associated 
challenges in the subsequent analysis of the LSS, and prepar- 
ing one for the analysis of the actual survey data. Further, 
the mock catalogues due to all the models and due to all the 
various biasing schemes were constrained to reproduce the 
two-point correlation function (Baugh 1996) and the power 
spectrum jBaugh fc Efstathioull993f) derived from the APM 
catalogue. The library of catalogues is impressive, both due 
to a wide coverage of the cosmological models which it of- 
fers, as well as in providing data due to a variety of biasing 
schemes, all consistent with the actual LSS as measured by 
the 2-point correlation function on scales between 1 to 10 
h~^Mpc. Further, the library provides 10 realizations for 
the rCDM model, which could be useful in assessing the ef- 
fect of cosmic variance on the statistics used to quantify the 
LSS. Due to these points, these catalogues are ideal to test 
various statistics which go beyond the usage of 2-point cor- 
relation function and which are likely to be sensitive to the 
higher order correlation functions. In fact, the sensitivity of 
all such statistics against various biasing schemes could also 
be tested. 



* We note however, that the SDSS redshifts are obtained using 
self-calibrated multi-band digital photometry, with the spectra 
obtained mainly in the Gunn-r band. Blanton et al.(2003) have 
obtained the luminosity function in all the bands using about 140 
000 galaxies from SDSS. 
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While 2dFGRS is complete and publically available and 
the SDSS team is making its releases annually, we take our 
present investigation in the spirit of "what to look for?^ when 
dealing with these datasets. Sheth et al.(2003) studied the 
geometry, topology and morphology of the dark matter dis- 
tributions in SCDM, rCDM and ACDM models simulated 
by the Virgo group. One of the motivations in carrying out 
this work is also to analyse as to what changes are brought in 
the geometry and topology of the excursion sets of rCDM 
and ACDM due to biasing^. By studying the constrained 
"galaxy distributions" in conic volumes of the mock cata- 
logues, we further hope to confront at least, some of the 
complications which one might anticipate in dealing with 
actual data. The present exercise will provide us with a the- 
oretical framework to compare the full SDSS data sets when 
these become available. 

O ur present study, in i ts scope, is similar to that carried 
out bv lCoUev et alj (|200(J) . However, unlike these workers, 
who studied the topology of the mock SDSS catalogues, our 
thrust is to test more complete set of statistics against the 
mock data. We shall of course be using more advanced and 
latest calculational tools for the purpose of our analysis. 

The remainder of this paper is organised as follows. Sec- 
tion 2 is devoted to a brief summary of the simulations which 
were used by CHWF to generate the mock catalogues. Here 
we also describe a subset of the mock catalogues which we 
use in our work. In Section 3, we elaborate upon the method 
used to extracting the volume limited samples from these 
mock catalogues for our further analysis. In Section 4 we 
study the percolation properties of the samples. Section 5 
briefly summarises the definition of MPs and the triangu- 
lation method used to model isodensity contours and for 
calculating the MFs for these contours. Section 6 outlines 
our method of analysis and presents our main results. We 
conclude in Section 7 by discussing the implications of this 
work and outlining its future scope. 



2 THE SDSS MOCK CATALOGUES 

2.1 Cosmological Models and their Simulations 

Cole, Hatton, Weinberg and Frenk (1998) (CHWF) prepared 
a comprehensive set of mock catalogues within a variety of 
theoretically motivated models of structure formation. Each 
of the models was completely specified in terms of a set of 
parameters defining the background cosmology and a set of 
parameters fixing the power spectrum of the initial density 
fiuctuations. 

The simulated models naturally fall into two cate- 
gories: (1) The COBE normalised models, where the am- 
plitude of the power spectrum is set by the amplitude of 
the fiuctuations measured in CMBR by COBE and extrap- 
olated to smaller scales using standard assumptions. (2) 
The "structure-normalised^ models wherein the amplitude 
of fluctuations is quantified in terms of erg, the rms density 
fiuctuations within spheres of radius 8 h~^Mpc. erg has been 
set from the abundance of rich clusters in our local Universe. 

The mock catalogues analysed in this paper are derived 

^ It should bo noted that the parameters for these models used 
by CHWF match with the models simulated by the Virgo group. 



from the structure-normalised models. Hence, we confine our 
present discussion only to this class of models. 

The description of the power spectrum is complete once 
the power spectrum index n and the Shape Parameter P 
are specified. In both the models of our intereset, n = 1 
(the Harrison- Zel'dovich power spectrum), and P has been 
chosen to be 0.25. ^ . This could be physically motivated ei- 
ther by having h = or by a change from the standard 
model of the present energy density in the relativistic par- 
ticles, e.g., due to a deca ying neutrino model proposed by 
iBond fc Efstathioul Jl99lD . 

The rCDM model is similar to standard CDM in hav- 
ing 0.0 = 1, but it has the same power spectrum as other 
low-density fiat models. For the sake of testing the discrim- 
inatory power of the MP-based statistics, we also choose 
to apply our methods to a realization of ACDM structure- 
normalised flat model with f2A=0.3. The ag in both the mod- 
els was_^xed__b^_th«_2resci^^ as provided 
bv IWhite. Efstathiou fc PrenkI lll993h . This corresponds to 
al = 0.55 and ag = 1.13'^. 

Starting with the complete specification of the model- 
parameters, the initial conditions of all the simulations were 
set according to glass-configuration and the simulations were 
evolved with a modified version of Hugh Couchman's Adap- 
tive PPM code llCouchmanlll99j) . The physical size of the 
box used was 346.5 h~^Mpc and the calculations were car- 
ried out on a 192^ grid. The simulations were evolved up to 
the present epoch. We refer the reader to CHWF for further 
details. 



2.2 Selecting the formation sites for galaxies 

The simulations described above provide the distribution of 
dark matter. In order to prescribe the preferred sites for 
the formation of galaxies within a given mass-distribution, 
CHWF employed a set of models which invoke a local, 
density-dependent bias. The rationale behind this exer- 
cise comes from t he ob servational evidence obtained by 
iPeacock fc Dodd3 il994ft . who found the distributions of 
Abell clusters, radio galaxies and optically selected galaxies 
to be biased relative to one another. The biasing schemes 
employed by CHWF were in terms of simple parametric 
functions because at the time, it was still not possible to de- 
termine the function that relates the probability of forming a 
galaxy to the properties of the mass density field. It should 
be noted however, that the situation has changed in the 
mean time due to development of elega nt analytic tools like 
the Halo Distribution Function (HDP) iBerlind et alJl2003l : 
iKravtsov et~ai]l2003l) and Conditional Lumino sity Function 
fCLFI iYang. Mo. van den Bosch fc Jindl2003l) which pro- 



^ We note that the current data due to WMAP constrain the 
value of r at around 0.21 which is quite close to the value adopted 
in models which we analyse in this paper. 

^ The current value of (Tg as provided by WMAP is (78=0.84. 
While our line of approach may remain the same, we hope to com- 
pare the observations with more realistic set of mock catalogues 
involving improved values and bounds on the above parameters. 
A detailed comparison with the observations is, hence, deferred 
to a future work. 
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vide us with a deeper understanding of issues relating to 
bias. 

CHWF ascribed the values of the parameters in the 
parametric, bias-prescribing functions by constraining the 
distributions of "galaxies" within each of their cosmological 
simulations to reproduce the two-point correlation fu nction 
on the scales between 1 to 10 h"^Mpc jBaughlllQQ^') . The 
ACDM and rCDM models which we investigate in this pa- 
per were populated with galaxies by following a selection 
probability function given by 

P{v) oc exp{aiy + l3u^^'^) if ^ 

oc exp(Qfi^) if ^ 0, (1) 

where the dimensionless variable v is defined to be i/(r) = 
5s{r)/as- Ss{r) is the density contrast of the initial den- 
sity field smoothed using a gaussian window on a scale of 3 
h~^Mpc. as is the r.m.s. mass-fluctuation of this field. For 
ACDM model, aA=2.55 and /3a=-17.75. The large negative 
value of /3a indicates that this model requires large anti-bias 
in the high density peaks of the initial density field in or- 
der to reproduce the observed 2-point correlation function. 
For the rCDM model, q:^=1.10 and /3r=-0.56. The above 
probability function was normalised so as to produce a total 
of 128^ galaxies within the cubic box of size 346.5h-iMpc, 
which corresponds to the observed average number density 
of galaxies with luminosity ^ L,/80. Further details about 
this and other biasing schemes can be found in CHWF. 



2.3 Preparing a mock catalogue 

Although CHWF have prepared mock catalogues for both 
2dFGRS and for SDSS, we shall confine our present discus- 
sion to the SDSS catalogues alone. 

When completed, the SDSS shall cover 3.11 steradian, 
i.e., about 1/4"^ of the sky. The survey geometry of the 
SDSS is defined in terms of an ellipse in the northern galactic 
hemisphere with the centre of the ellipse at R.A.= 12*^^20™, 
5=32.8°, which is in close proximity of the northern galactic 
pole. The semi-minor axis of the ellipse runs for 110° along a 
line of constant R.A., whereas the semi- major axis subtends 
an arc of 130°. CHWF choose the location of an observer 
inside the cube, and mimic the survey geometry around the 
observer's position. The selected galaxies are taken to fall in- 
side the survey geometry prescribed above. The simulation 
cube is replicated periodically in all the directions to achieve 
a depth corresponding to z=0.5. For the sake of ideal com- 
parison, the observer's position in all the mock catalogues 
is chosen to be the same. 

The SDSS operates on its own multi-band digital pho- 
tometry and primary selection of galaxies is made in Gunn-r 
band. However, for the sake of simplicity, CHWF prepare the 
mock SDSS catalogues using a hj band luminosity function 
with a Schecter form. By doing so, they are able to reproduce 
the radial selection function and the expected number count 

* On observational side, lYan. Madgwick fc Whitd j20o3) 
have constrained the halo m odel using 2dFGRS whereas, 
IMagliocchetti fc Porcianil i2003l) study the halo distribution of 
2dFGRS galaxies. At least in within a given cosmological model, 
hence, a robust study of bias should now be possible. 



for the SDSS ^ . The magnitude limit employed is bj <18.9, 
so as to approximately match the expected SDSS count of 
900 000 galaxies in the survey area. The Schecter luminos- 
ity function, which gives the comoving number density of 
galaxies with luminosity between L and L-l-dL, is given by 

^(L)dL = |l exp(-i/L.)dL, (2) 

where CHWF employ 0. = 1.4x 10"^h^Mpc"^ a, = -0.97 
and M^j-Slogioh = -19.5. 

The luminosity function and the magnitude limit com- 
pletely specify the radial selection function, which prescribes 
the number of observable galaxies as a function of radial 
distance r from the the observer (assuming isotropy and no 
evolution in the luminosity function over the redshift-range 
of interest; z<0.5). This information is used to select an ap- 
propriate number of galaxies in a shell of thickness dr at a 
distance r. These galaxies are further randomly assigned lu- 
minosities consistent with the Schecter form of the luminos- 
ity function. Having ascribed the luminosity, CHWF also tag 
every galaxy with the maximum redshift Zmax up to which it 
would be detectable in the survey. As it turns out, the pro- 
vision of this number makes it easier for the user to generate 
a volume limited sample out of a given mock catalogue. 



3 DATA REDUCTION: EXTRACTING THE 
MOCK SAMPLES 

3.1 Preparing Volume Limited Samples 

With this section, we return to the subject of our present 
investigation. Our interest lies in utilising the geometric and 
topological characteristics of the excursion sets of a cosmic 
density field to test whether we can compare and distinguish 
between the rival models of structure formation. For this 
purpose, we work with 10 realizations based on the mock 
SDSS catalogue due to rCDM model and compare it with a 
realization of the ACDM model. 

Since the excursion sets of a density field are extended 
objects, one should ensure that no user-specific bias is intro- 
duced while studying them. However, precisely such bias is 
encountered in dealing with flux-limited samples. Such sam- 
ples are characterised by a selection function which falls off 
upon moving radially outwards. This is because the galaxies 
become fainter with increasing distance, and the cut-off lu- 
minosity Lmin,(z) corresponding to the apparent magnitude 
limit of the survey increases with redshift, leading to a drop 
in the number of detectable galaxies. 

To illustrate our point consider the case of an extended 
filamentary supercluster (which could contain up to 10'*"^ 
galaxies) and which is extending radially outward from the 
location of the observer. In a ffux limited sample the selec- 
tion function will ensure that galaxies of a given brightness 
belonging to the rear end of the supercluster are system- 
atically suppressed in number compared to similar galaxies 
closer by. Thus the size and shape of the supercluster will 
be grossly distorted in a ffux limited sample. It is clear from 

^ The bj band luminosity of a galaxy is related to the SDSS g' 
and r' bands: bj = g'-|-0.155-f 0.152(g'-r'). 
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here that the physical definition of a supercluster here de- 
pends on the position and orientation of the supercluster 
relative to the observer. In principle, no such dependence 
should be allowed in defining the supercluster. 

To avoid such problems and to ensure that our analysis 
is bias-free, we shall first construct volume limited subsam- 
ples from the fiux limited mock surveys and then apply our 
methods to derive morphological and geometrical properties 
of these volume limited subsamples. 

This is achieved by discarding all such galaxies at a 
given distance r, which fall below the detection limit at the 
farthest end of the sample. In turn, only such galaxies are 
selected from the flux-limited sample, which would be visible 
throughout the volume of the sample being analysed. This 
is equivalent to rendering the selection function spatially 
invariant. 

The selection function can be written as 

/oo 
(t>{L)dL, (3) 

where di is the luminosity distance corresponding to ACDM 
and rCDM, where appropriate. The quantity L,„i,i(z) cor- 
responds to the apparent magnitude limit of the survey and 
is an increasing function of z. Figure Q shows the selection 
function S(z) along with the survey-volume in (h~^Mpc)"^ 
for both ACDM and rCDM models. The selection function 
S(z) stands for the number of galaxies included in the sur- 
vey within a spherical shell of radius dL(z) and an infinites- 
imal thickness Ad_L. In both models S(z) peaks at a redshift 
~0.11 and falls off on either side. The fall-off of S(z) on the 
left side is a volume effect, i.e., though the number density 
of detectable galaxies is higher (than at z~0.11), the avail- 
able survey-volume is small. For z>0.11, the number density 
of galaxies is lower than at z=0.11, even though the physi- 
cal volume is large. From the point of view of analysis, one 
would prefer as large a volume as possible, and a constant 
number density of galaxies to pick out from all parts of this 
sample- volume. 

Choosing the depth corresponding to z>0.11 increases 
the survey volume at the cost of the number density. The 
mean inter-galactic separation A(~ n'^^^) in such sample(s) 
will be larger. Hence, to obtain a continuous density field, 
we need to smooth the galaxy-distribution more. For z<0.11, 
the inter-galactic separation is reduced; we need to smooth 
less, but this happens at the cost of the survey volume. 
Hence, a situation may be envisaged wherein the structures 
which we study do not represent the large scale distribution 
of matter in a fair manner. Thus, the peak of S(z) signi- 
fies a compromise between two competing factors, the num- 
ber density and the survey volume. By working with the 
volume-limited sample corresponding to the peak of S(z), 
we optimise on the inter-galactic separation (and hence, the 
required smoothing) and on the sampling of the large scale 
matter-distribution (which is proportional to the available 
survey volume.). 

From the above analysis we expect to achieve a fair 
sample of LSS if we work with a volume-limited sample 
of depth z~0.11. Creating a volume-limited sample out of 
the mock catalogues of CHWF is relatively simpler because, 
in addition to supplying the luminosity and the redshift of 
the galaxy, CHWF also provide the maximum redshift Zmax 
upto which the galaxy will be visible. In order to prepare a 



volume- limited sample with a given depth Rmax, we select 
aZZ the galaxies which are nearer to us than Rmax, and which 
have Zmax ^ Zmax, the limiting redshift of the sample. We 
study the number of galaxies included in the sample as a 
function of limiting redshift Zmax- While we observe this 
number to peak at z~0.11 in both the models (Figure Ql, 
we prefer to work with rCDM samples of depth '^ma'x~ 

0.117 

(Ng ~2. 19x10*). Ng, the number of galaxies in rCDM, is 
comparable to the number of galaxies included in ACDM 
catalogue if Z^^x— O-H- The solid and dashed vertical lines 
in both panels of Figure^refer to these redshifts for rCDM 
and ACDM models respectively. Our choice of limiting red- 
shifts is motivated by our intention of working with same 
physical volumes for both the models. The horizontal line in 
the right panel of this figure shows that the survey volumes 
in both models are comparable. We finally have 10 rCDM 
samples and a ACDM sample, all having practically iden- 
tical volume and a similar number of galaxies. This means 
that the mean inter-galactic separation is similar in all 11 
samples and we can employ the same smoothing scale in all 
cases. This is essential in the wake of assessing the regularity 
behaviour of MFs and to compare the two rival models at 
hand. 

The mock catalogues provide us with redshift-space po- 
sitions of the galaxies in terms of their Cartezian coordi- 
nates {x,y,z). The Z— axis of the coordinate system points 
towards the centre of the SDSS survey. The longitude I 
{—tan~^{y/x)) is measured relative to the longer axis of the 
SDSS ellipse, whereas the latitude 6=90° is the centre of 
the SDSS survey. Knowing {x, y, z), it is straightforward 
to obtain the Cartezian coordinates (X, Y, Z) of the galax- 
ies on a distance scale. Throughout we work with redshift- 
space coordinates of the galaxies , and do not attempt to 
correct for the redshift space distortions. This decision is 
partly influenced by the analysis of CoUey et al.(2000), who 
found the global genus-curve of the mock SDSS catalogue 
to be only mildly sensitive to redshift-distortions. Since the 
redshift space distortions are directly sensitive to Qo, we 
do anticipate definitive signatures of redshift space distor- 
tions on the sizes of the structures. In particular the fin- 
gers of God effect needs to be quantified using MFs and 
Shapefinders. Such an investigation is beyond the scope of 
the present paper and will be taken up in a future work. 
The typical box enclosing all the samples is characterised by 
)~(--322, — 290.5,0.0) and (X^ax, Ymax, 
Zmaa;)— (322, 290.5, 357), where all the coordinates are mea- 
sured in h~^Mpc. Figure |5| shows the central slices of the 
ACDM and one of the realizations of rCDM mock cata- 
logues, with a similar number of galaxies. These slices refer 
to the same set of random numbers in the simulations. 



3.2 Generating Density Fields 

Our analysis will be carried out on a continuous density 
field sampled on a grid with uniform resolution. Thus, hav- 

CHWF assume that the "galaxies" in their mock catalogues 
share the bulk flow with the matter distribution, and do not at- 
tempt to model the velocity field of galaxies in any more detail. 
Their decision was at the time supported by the Virgo simulation 
studies by Jenkins et al.(1998). 
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Figure 1. Left panel: The selection function is shown for both the models. S(z) peaks around z=0.11 in both the cases and falls off on 
either side. The solid vertical line refers to the limiting redshift (z=0.117) for the tCDM volume limited samples. The dashed vertical 
line signifies the limiting redshift (z=0.11) for the ACDM model. All the samples include ~2. 19x10^ galaxies, which corresponds to a 
mean inter-galactic separation of ~6 h~^Mpc. Right panel: The survey- volume is plotted as a function of the limiting redshift z. We 
notice that the limiting redshifts (solid and dashed vertical lines) refer to similar values for the survey-volume. 



ing extracted the volume limited samples, we next present 
our method of generating continuous density fields for these 
samples. 

As noted in the last subsection, the size of the box en- 
closing a typical sample is (AX, AY, AZ)~(644, 581, 357) 
h"^Mpc. We fit a grid of size 184x166x102 onto this box. 
The size of the cells, i.e., the resolution, turns out to be 
3.5 h~^Mpc. About 63 per cent of the cells are found to be 
inside the survey boundaries. We first carry out a Cloud- 
in-Cell (CIC) smoothing by attaching suitable weights to 
each of the eight vertices of the cube enclosing the galaxy. 
This procedure conserves mass. About 25 per cent of the 
cells inside the survey volume acquire nonzero density after 
CIC smoothing. As we noted earlier, the mean inter-galactic 
separation is A ~ n^^''^ ~6 h~^Mpc ~ 1.71x£g, where Ig 
is the resolution of the grid. Thus we do expect the CIC- 
smoothing to leave behind an appreciable fraction of cells 
with zero density. Ideally, a density field should be continu- 
ous, and should be defined over the entire sample volume. To 
achieve this, we need to smooth the CIC-smoothed density 
field further. We employ a gaussian window 

W^4r) = ^exp(-^) (4) 

for this purpose. A gaussian window, being isotropic, effec- 
tively smooths the matter distribution within As h~^Mpc 
of the cell. Since the purpose of our present investigation 
is to quantify the morphology of coherent structures, the 
smoothing scale As should not be very large. It should be 
just enough to connect the neighbouring structures. While 



iGott et alj lll989tl prescribe using As ^2.5^9 for such appli- 
cations we prefer to work with As=6 h~^Mpc, which is 
comparable to the correlation length scale in the system, 
and reasonably smaller than 2.5^^. As we shall demonstrate 
below, the behaviour of MFs is remarkably regular even at 
this smoothing scale. The lower panel of Figure |5| shows 
the adopted smoothing scale in the overall distribution to 
develop a feel for the effective smoothing done with the 
samples. It is clear from here that the large voids should 
be practically devoid of matter even after smoothing and 
should be accessible for morphological analysis at the low- 
est, sub-zero 5- values, whereas the nearby clusters and su- 
perclusters should connect up at a suitable density thresh- 
old (corresponding to (5~l-2, see below). The morphology of 
these superstructures should be accessible at or near the on- 
set of percolation. In this paper, we shall concentrate on the 
shapes, sizes and geometry of such superclusters and shall 
quantify these notions in more detail in Section 6. 

Prior to smoothing by a gaussian (or say, a top-hat) win- 
dow, which requires utilising FFT routines, a suitable mask 
for the cells outside the survey boundaries has to be as- 
sumed. It is also to be noted that smoothing leads to a spill- 
over of mass from the boundaries of the survey to the outer 



An isotropic smoothing with sufficiently large As (say, 
^10h~^Mpc) will reduce our discriminatory power by diluting 
the true structures to a great extent. On the other hand, if one 
intends to compare the MFs with their theoretical predictions 
(Matsubara 2003), one should employ large scale smoothing in 
order to reach a quasi-linear regime. 
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Figure 2. Shown here are the central slices of the ACDM mock 
SDSS catalogue and of one of the realizations of its rCDM coun- 
terpart. Both the catalogues refer to the same set of random 
numbers. The coordinates refer to the redshift space positions. 
The fingers-of-God effect is evident from these slices. The galaxy- 
distribution is constrained to match the observed 2-point correla- 
tion function, and the two slices qualitatively look quite similar. 
The voids in ACDM are however noted to be much cleaner, and 
the fingers-of-God effect stronger. The galaxy distributions are 
smoothed with a gaussian window of size 6 h~^Mpc. In the left 
corner of the bottom panel we show the window-size to give a feel 
for the effective smoothing done around every galaxy. 

peripheries. Thus, the mass inside the survey is no longer 
equal to the mass contained in all galaxie s. To conserve mass, 
we em ploy the smoothing prescription of lMelott fc DominikI 
il993l) : we prepare two samples; the first containing the CIC- 
smoothed density field, with cells outside the survey volume 
set to zero density, and the second reference grid wherein the 
cells inside the survey volume are all set to unity, while ones 
outside the survey boundaries are set to zero. Both density 
fields are smoothed with the gaussian window (Eq|lJ at As 
— 6h~^Mpc The reference grid helps us correct for the 
effect of having a smoothing window at the boundaries of 
the survey. We divide the density values of the first grid, 
cell-by-cell, by the entries for the density values in the ref- 
erence grid. We find that this corrects for the spill-over of 
mass at the boundaries, and we indeed find mass inside the 
survey volume to conserve to a very good accuracy. 

Having elaborated on our method to prepare volume 

We use the Fastest Fourier Transform in West (FFTW) rou- 
tines for the fourier transforms carried out on our grid. Since the 
grid-dimensions here are not powers of 2, standard fourier trans- 
form routines cannot be used for this purpose. 



limited samples and to recover the underlying density fields, 
we now elaborate upon our method of analysis and finally 
present our results. Before embarking on quntifying the mor- 
phology of density fields, we shall study the percolation 
properties of the samples. This is the subject of the next 
section. 



4 PERCOLATION OF THE LARGE SCALE 
STRUCTURE 

Percolation gives us useful insights into the connectivity 
of the density fields. Earlier workers (for example, Sahni, 
Sathyaprakash & Shandarin 1997) have found noticeable 
correlations between the topological behaviour of density 
fields and their percolation properties. Density fields show- 
ing bubble topology, meatball topology or sponge/network 
topology have distinct features in the way such fields per- 
colate. The 3-D idealized gaussian random fields (here- 
after, GRFs), which often serve as a useful benchmark 
in the study of large scale structure, percolate at 16 per 
cent of the overdense/under dense volume filling fraction 
JShandarin fc Zeldovichlll989t) . and those defined on a grid 
percolate at ~ 30 per cent volume frac tion for a variety of 
power spectra jYess fc Shandariiilll996r) . Overdense regions 
undergoing gravitational clustering percolate at lower values 
of the filling fraction as compared to GRP. As a result, the 
difference between the thresholds of percolation (measured 
in terms of the volume filling fraction) serves as a useful 
diagnostic with which to quantify the degree of clustering. 
Such fields show network topology and are characterised by 
"clusters" having larger length-scales of coherence. Various 
workers (e.g., Sheth et al., 2003) also find that the dom- 
inant morphology of the superclusters (defined for a field 
when the percolation sets in) is filamentary. Thus, study- 
ing the percolation properties of our present catalogues can 
help us identify a threshold of density at which to look for 
realistic structures in our survey volumes. 

In the ensuing analysis, our results shall refer to 11 den- 
sity fields (10 mock catalogues due rCDM and one due to 
ACDM) scanned at 50 thresholds (or levels) of density. This 
sampling is not uniform. Our preliminary analysis shows 
that overdense regions percolate when the overdense vol- 
ume fraction ^ 10 per cent of the survey volume. On the 
other hand the underdense regions or voids percolate when 
the overdense volume fraction is ~70 per cent. Thus, from 
a morphological perspective, the volume fractions 5C30 per 
cent and ^70 per cent are interesting for overdense and 
underdense regions respectively. In the intermediate range, 
both the overdense and underdense regions percolate, and 
most of the volume is occupied by the largest cluster and the 
largest void respectively. So practically very little is gained 
by studying the structures in the range FFv £(0.3,0.7), 
where FFy, the overdense volume fraction is given by 

= ^ J ©(P - PTH)d\. (5) 

Here 0{x) is the Heaviside Theta function. 

e{x) = 1x^0 

= X <0. (6) 

FFv measures the volume-fraction in regions which sat- 
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Figure 3. The density contrast and tlie mass filling fraction FFm are plotted as functions of FFv, tlie volume filling fraction which is 
treated as an independent parameter. The solid line stands for the tCDM model and is obtained by averaging over 10 realizations. The 
dashed line stands for the ACDM model. Please refer to the text for further comments. 



isfy the 'cluster' criterion pciustcr Pth at a given density 
threshold Pth- Vt is the total survey volume. In the fol- 
lowing, we use FFv as one of the parameters with which 
to label density contours. Finally, the ranges FFv ^0.3 and 
FFv each are attributed 23 levels of density, each equi- 
spaced in FFv. The intermediate range FFv €(0.3,0.7) is 
sampled at four levels all equi-spaced with AFF=0.1. In Fig- 
ure |3 we study the density contrast S and the mass filling 
fraction FFm with respect to FFv^'^. The first noticeable 
effect which the bi asing brings in is, that unlike the dark 
matter simulations JSheth et alJl2003D . FFm no longer dis- 
tinguishes the two models. At all the chosen levels, the over- 
dense mass fraction, i.e., the fraction of galaxies enclosed by 

Throughout our discussion, we take the mass enclosed within 
a supercluster to be equivalent to the number of galaxies enclosed, 
since the two differ from each other only by the mass of the parti- 
cle in the simulations, which is a constant. Though this is different 
for the two models, the mass filling fraction which usually enters 
our discussion, is unaffected by this. 



a given normalised volume FFy, is similar in ACDM and 
rCDM. The mass-parametrisation which was noted to be so 
effective in distinguishing dark matter distributions in vari- 
ous models, cannot perform any better than FFv- We shall 
return to this issue in Section 6. We further note that, the 
density contrast S follows similar pattern in both the models. 

We employ a friends-of-friends (FOF) based cluster- 
finding code to study the number of clusters and the over- 
dense volume fraction in the largest cluster (FFlc) at all the 
50 density- levels. FFlc stands for the volume of the largest 
cluster normalised by the total overdense volume. 

For cubic samples, the percolation is said to have set in 
when the largest cluster first spans across the full length of 
the box (by periodicity arguments, this is equivalent to hav- 
ing the supercluster spanning across the infinite space). In 
the case of redshift surveys wherein the matter-distribution 
generally falls within a cone, we cannot utilise such a con- 
vention, for the survey boundaries diverge as we go radi- 
ally outwards and a supercluster spanning the survey vol- 
ume could as well be touching the two boundaries nearer 
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Figure 4. For samples with arbitrary survey geometry, the growth-rate of the volume of the largest cluster is important for quantifying 
the onset of percolation. Here this quantity is plotted w.r.t. FFy. The merger activity in tCDM is pronounced for FFv6[0.06,0.1]. We 
take FFv ~0.08 to be the percolation threshold for this model. Similarly, in case of ACDM, the largest cluster growth is highest at 
FFv '-0.08. Hence, we take FFp^^.;^ — 0.08. Where needed, we explore the whole range of density levels <5G[<5pcrc,<5max]. 



to the observer, without visiting other parts of the survey 
volume. Clearly, this does not mean the system has perco- 
lated. Thus, for conic survey volumes, spanning the survey 
volume is not a robust concept to signify the onset of percola- 
tion. Instead, Klvuin & Shandarin (1993) motivate utilising 
the merger property of clusters to define the onset of per- 
colation. Following this convention, we measure the rate of 
growth of the volume of the largest cluster. The highest rate 
of growth signifies pronounced merger of clusters, leading to 
a single, larger cluster. The threshold of density at which 
the growth-rate is highest is defined to be the percolation 
threshold. Figure |1| shows the rate of growth of the largest 
cluster as a function of FFy. For rCDM, we use the data 
after averaging over all 10 mock realizations, and find that 
the merger rate is highest in the range FFv€[0.06,0.1]. For 
our purpose, we take FFy of 8 per cent to be the percolation 
threshold for rCDM model. We find that the largest clus- 
ter of ACDM undergoes a pronounced merger at FFv~0.8, 
bringing in 30 per cent of the overdense volume fraction 
into the largest cluster, (see Figure |^ . We have only one 
realization for the ACDM model, so a statistically robust 
statement cannot be made, but for all purposes, we treat 
FFv= FFP'"'''=0.08 as the percolation threshold for ACDM. 



Where needed, we explore the entire range of FFv ^0.4 for 



various statistics. To summarise, FF\ 



iFF^ 



=0.08. 



Finally, Figure |S] shows the percolation properties of 
our samples. We study the cluster abundance (lower pan- 
els) and the filling fraction in the largest cluster FFlc (the 
upper panels) as functions of three different parameters, 5, 
FFv and FFm. The parameter values corresponding to per- 
colation are marked with a dotted line for both the models. 
For rCDM, we again show the curves after averaging over 
the 10 available realizations. We note that the clusters are 
slightly over-abundant in ACDM compared to rCDM. As- 
suming a similar degree of scatter in the cluster-abundance, 
we conclude that the Number of Clusters statistic is useful 
to discriminate the two models. The quantities FFv^'^'^'^ are 
quite closer to each other for both models as we noted ear- 
lier (middle, upper panel). From here, we conclude that the 
spatial connectedness, as measured in terms of percolation 
of the volume of the LSS alone is not sufficient to distinguish 
between the two models. The top-right panel does indicate 
that in rCDM, 30 per cent of the total mass is in "clusters" 
when percolation sets in, and 32 per cent of this mass (which 
amounts to ~22000 galaxies\) is in the percolating superclus- 
ter. Looking at the similar volumes that these systems oc- 
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Figure 5. We study here the cluster abundance and FFlCi the fraction of overdense volume in the largest cluster, with three different 
parameters, S, FFy and FFm- The rCDM curves are averaged over 10 realizations. The percolation of the volume of the largest cluster is 
remarkably similar in both the models. The percolating supercluster occurs at (5pci-c— 1.7 in both ACDM and rCDM. The dotted vertical 
lines in all the plots signify the onset of percolation. We also note from the lower panels, that the number of clusters in ACDM are larger 
than in rCDM. This could be a useful means to distinguish these two models. 



cupy at percolation (although the ACDM percolation sets in 
a little earlier than in case of rCDM) , clearly other geomet- 
ric and topological descriptors like Minkowski Functionals 
would be more insightful to assess the properties of the su- 
perclusters. In passing, we note two important aspects of the 
models: (1) The upper left panel shows that 5 is distributed 
over the same range of values for both the models. (2) For 



both the models, Spac (z[l,2j, which coincides with the con- 
ventional definition of superclusters. Thus, this range in S is 
in general useful for studying the network of superclusters. 
We shall make use of this fact in our subsequent analysis. 
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5 ESTIMATING MINKOWSKI FUNCTIONALS 

The work of iGott. Melott fc DickinsorJ dlQSd) led to a re- 
vival of interest in the paradigm of studying and quan- 
tifying the excursion sets of the cosmic density fields; 
an idea which was successfully propounded earlier in the 
topological context by Doroshkevich (1970). These early 
efforts however, were devoted to studying the topology 
and connectivity of isoden sity contours of the large scale 
structure (hereafter, LSSl llMelottI Il990l) . It was soon re- 
alised that topology alone may not be sufficient to de- 
scribe LSS, and that it should be complemented with ge- 
ometrical information pertaining to excursion sets. This 
goal was met in the form of the Minkowski Func - 
tionals (hereafter, MFsl dMecke. Buchert fc Wagner! I199I . 
A set of grid-based methods (e.g, the Crofton's for- 
mulae, the Koenderink Invariants) and a method to 
calculate MFs for poin t processes (the Bo olean grain 
model) were proposed l|Mecke. Buchert fc W agner 1994; 
ISchmalzing fc Buchertll997tl to address the issue of estimat- 
ing the MFs. (In fact, the CONTOUR3 D code to calculate 
the topology of LSS bv lWeinberd (^8^ is also grid-based.) 

Sheth, Sahni, Shandarin fc Sathyaprakash (2003) pro- 
posed a novel method to calculate MFs by following a math- 
ematically precise prescription with which to construct iso- 
density contours while quantifying them, as opposed to grid- 
approximations to contours implicit in earlier treatments. 
Their method accurately modelled isodensity contours in ad- 
dition to performing an online calculation of MFs for these 
contours. 

In 3— D, the MFs are four in number. The first three of 
these are geometric in nature, and change when we deform 
the surface, whereas the fourth is a topological invariant and 
tells us how many handles does a surface have in access of 
the holes which it completely encloses. The four MFs are 



• [1] Volume V enclosed by the surface, 

• [2] Area S of the surface, 

• [3] Integrated mean curvature C of the surface. 




• [4] Integrated gaussian curvature, or the Euler charac- 
teristic. 




A related measure of topology is the genus G — 1 — x/2, 
which is what we use in all our analysis. 

We shall evaluate MFs for mock catalogues by triangu- 
lating isodensity surfaces using SURFGEN. SURFGEN is 
a fusion of a robust surface modelling scheme and an effi- 
cient algorithm to calculate the geometric quantities. The 
algorithm to calculate MFs on a triangulated surface and 
the code SURFG EN have been discussed extensively by 
ISheth et alJ l)2003fl and we refer the reader to that paper for 
further details. On a dual processor DEC ALPHA machine, 
SURFGEN calculates the MFs for all the clusters found at 
about 100 levels of density on a typical 128"^ grid in within a 
few hours' time. This amounts to probing the MFs for about 
100 000 clusters of varying sizes. In fact, in our present ap- 
plication, which involves a grid which is 1.5 times larger than 
this, the calculations for all the 11 samples could be accom- 



plished in within 2 days' time. This highlights the scope of 
using SURFGEN in more ambitious applications involving 
larger number of realizations studied at higher resolution. 
Our method of analysis and the results are described in the 
next section. 



6 GEOMETRY, TOPOLOGY AND 

MORPHOLOGY OF MOCK CATALOGUES 

The mock SDSS catalogues being studied here, by construc- 
tion, have the same two-point correlation function in the 
range 1 to 10 h"^Mpc. Since Minkowski Functionals are 
shown to depend on the entire hierarchy of correlation func- 
tions, they may be expected to differentiate between models 
as similar as rCDM and ACDM. This will indeed be the 
case, as we shall show below. 

One of the fascinating aspects of the LSS is related to its 
visual impression. The LSS as revealed in various dark mat- 
ter simulations and redshift survey slices reveals a cosmic 
web of interconnected filaments running across the sample, 
separated by large, almost empty regions, called voids (see 
Figure 121 for an illustration). Our study will focus both on 
global MFs which characterise the full excursion sets of the 
cosmic density field and partial MFs which, together with 
Shapefinders, help us study the shapes and sizes of the struc- 
tural elements (superclusters) of the cosmic web 

We begin by first studying the global Minkowski func- 
tionals. Next we focus on how the largest cluster evolves 
with the volume fraction. The last two subsections deal with 
the morphology of individual superclusters in rCDM and 
ACDM models. 

6.1 Global Minkowski Functionals 

As mentioned earlier, we sample the 11 density fields (1 
ACDM + 10 rCDM) at 50 levels of density. The criteria 
adopted for their choice has been described earlier (see Sec- 
tion 4). The levels uniformly correspond to the same set 
of volume fractions FFv, for all the samples, which makes 
comparison between models easier. 

We find clusters using the friends-of-friends (FOF) al- 
gorithm based Clusterfinder code at every level of density. 
Next we run SURFGEN on all the clusters and estimate 
their MFs by modelling their surfaces. These are the partial 
MFs and are useful for studying the morphology of individ- 
ual clusters'^^. We use the additivity property of the MFs 
and evaluate global MFs by summing over the contribution 
of partial MFs due to all the clusters at a given level of den- 
sity. For rCDM, we use the measurements from available 10 
realizations to produce average global MFs. Figure |H| shows 
the dependence of the global MFs on the volume fraction 
FFv . The solid curves represent the mean global MFs due to 
rCDM model and the error-bars stand for the la standard 
deviation. The dashed curves stand for the ACDM model 

^■^ We do not however attempt to study the morphology of voids, 
which is an equally important constituent of the cosmic web. 

It should be noted however that the "clusters" referred to in 
the present discussion do not coincide with the galaxy-clusters. 
Our clusters are connected, overdense regions and at appropriate 
level(s) of density may coincide with the superclusters of galaxies. 
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Figure 6. Global Minkowski functionals are evaluated at 50 density levels at a common set of volume filling fractions and are studied 
here w.r.t. FFv- We show the values of global MFs as normalised to the volume of [100 h~^Mpc]''. The global MPs for rCDM model 
are averaged over the available 10 realizations (solid lines) and the errorbars represent la deviation. Assuming the same level of accuracy 
for ACDM, we may conclude that the global MFs with volume parametrisation can indeed clearly distinguish tCDM from ACDM. 
Considering the fact that all the samples show the same two— point correlation function, the above difference can be attributed to the 
higher order correlation functions. MFs are able to respond to such a difference, and hence can be taken to be robust measures for 
quantifying the cosmic web. For further discussion, please refer to the text. 



. First, we note that the global MFs of the rCDM model 



^° It should be noted here that our samples do not exhibit 
periodic boundary conditions. Hence, as explained in detail by 
ISheth et all |20o3), our measurements cannot be compared with 
the theoretical predictions straightaway, which are generally for 
periodic samples. In fact, the present measurements refer to the 
overdense excursion sets, and should be complemented by those 
for voids, to finally get results for an equivalent periodic sam- 
ple. See Shandarin, Sheth & Sahni (2003) for further details. The 
boundary effects in the present measurements are evident at large 



are extremely well constrained. This shows the remarkable 
accuracy with which SURFGEN will enable us to measure 
the MFs of the LSS due to large datasets like SDSS and 
2dFGRS. We further note that the two models can be dis- 
tinguished from each other with sufflcient confidence using 
the global area, mean curvature and genus measurements at 
FFv^O.5 whereas the global volume scales linearly with 

FFv, where the area and curvature saturate to nonzero values 
referring to an all-pervading, cone-filling "cluster" . 

This is to be expected, because 50 per cent volume fraction 
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FFv, as expected. Distinctly difFerent trends of MPs in two 
models are due to the contribution of the higher order corre- 
lation functions, and do indeed reflect the manner in which 
the galaxies cluster in ACDM as against rCDM. This can 
be dubbed to be the success of MFs as robust statistics ideal 
for quantifying the cosmic web. 

From the dynamical poi nt of vi e w, the s e results carry 
some surprises. As noted bv lMelot3 il990l) : ISpringel et alJ 
l|l998h . the amplitude of the genus-curve drops as the 
N— body system develops phase correlations. Given two den- 
sity fields, the system with larger genus-amplitude shows 
many more tunnels/voids which are, therefore, smaller in 
size. As time progresses, the voids are expected to expand 
and merge, leading to a drop in the genus-amplitude, while 
the phase correlations continue to grow. With this simple 
model, one could correlate the amount of clustering with 
the relative smallness of the amplitude of genus, and there- 
fore, of the area a nd the mean curvature. As illustrated by 
ISheth et all ll2003ll . the dark matter distribution of rCDM 
model due to Virgo group consistently shows considerably 
larger amplitudes for the MFs compared to the ACDM 
model whereas, we find the reverse trend in the MFs 
of "the galaxy distributions" due to the same two mod- 
els. This is a sufficiently robust result because, the rCDM 
curves are averaged over 10 realizations. We could inter- 
pret this as an instance of the phase-mismatch between the 
galaxy-distribution and the underlying dark matter distri- 
bution and could attribute it to the biasing prescription in- 
voked. This is supported by the fact that CHWF require 
large anti-bias in high density regions of the ACDM model. 

To conclude, the study of the global MFs reveals 
from here that, local, density-dependent bias could lead to 
an apparent phase-mismatch or different phase-correlations 
among the dark matter and the galaxies. In passing, we also 
note that ACDM preserves its unique bubble shtft which is 
pronounced compared to rCDM (lower right panel of Figure 

Figure|7|shows the global MFs plotted against the over- 
dense mass-fraction FFm. The mass- volume relation is re- 
markably similar among both the models. This again is an 
effect of biasing. In addition to the fact that the models 
can be distinguished, we also note that the excursion sets 
enclosing the same fraction of mass due to the two models 
show appreciably different geometry and topology. Thus, as- 
sessing the excursion sets enclosing the same mass-fraction 
could be a useful means of analysing the realistic data and 
comparing them with the mock data. 



6.2 Minkowski Functionals for the Largest Cluster 

As we scan the density fields by lowering the density thresh- 
old, the largest cluster of the system grows in size due to 
merger of smaller clusters. The manner of evolution of the 
largest cluster is connected with the percolation properties 



refers to ~85 per cent of the mass fraction (Figure |2J. Thus, the 
excursion sets at FFy^O.S already contain most of the mass, and 
the necessary geometrical and topological information is conveyed 
by studying them. 

However, rvirgo=0.21, which is different from r=0.25 adopted 
in the simulations studied here. So long as F is the same for both 
the models, the conclusions are expected to remain unaltered. 



of the system. We noted earlier (Section 4) that the largest 
cluster spans most of the volume by the time the volume 
fraction reaches FFv~0.3. In fact, most of the features in 
the global MFs beyond FFv=0.3 pertain to the evolution of 
the largest cluster alone. Hence, in this section, we study the 
evolution of the MFs of the largest cluster for FFy^O.S. We 
find that, of the four MFs, the integrated mean curvature 
(hereafter, IMC) and the genus are the best discriminants. 
Further, it is instructive to see the behaviour of the largest 
cluster with FFm— parametrisation. In Figure |H| we show 
the IMC and genus for the largest cluster as functions of 
FFv (left panels) and FFm (right panels). We notice in the 
bottom-left panel of this figure that the genus of the largest 
cluster remains vanishingly small until the onset of perco- 
lation in both the models. However, once percolation sets 
in, the genus of the largest cluster grows much more rapidly 
in ACDM than in rCDM. The larger genus is indicative of 
greater number of tunnels (more porosity) which the ACDM 
structure exhibits relative to its rCDM counterpart. While 
the genus in case of rCDM seems extremely well-constrained 
(small scatter), the fluctuations in IMC, at least at smaller 
FFy— values (<0.15), are reasonably large. This could be in- 
ferred as fluctuation in the sizes of the tunnels forming due 
to the merger of the smaller clusters, while their total num- 
ber conserves for a given FFv. The IMC and genus of the 
largest cluster in the range FFv€[0.05,0.3] could be some of 
the best discriminants of the two models. 



6.3 Morphology of the Superclusters 

This and the following subsection are devoted to a critical 
study of the morphology of the ensemble of superclusters. 
Based on a statistical study of many realizations of rCDM, 
we are able to establish it for the first time that the super- 
cluster shapes and their sizes are very well constrained mea- 
sures, suitable as ideal diagnostics of the large scale struc- 
ture of the Universe and are sensitive to the cosmological 
parameters of the model(s). 

Since we are dealing with density fields, our superclus- 
ters are defined as connected, overdense objects. In Section 
4, we studied the percolation properties of the density fields 
and found that the percolation takes place at a threshold of 
density corresponding to 5=1.7 in both ACDM and rCDM. 
Further, Jmax — 3.5. Thus at the highest threshold of den- 
sity, we are still within the weakly nonlinear regime, and may 
find sufficiently large connected objects which might act as 
dense progenitors of the percolating superclusters. With this 
anticipation, we find it advisable to study the morphology 
of the ensemble of superclusters as a function of the density 
contrast over the entire range 5G[0,5max]. In so doing, we are 
also able to assess how the morphology of the LSS changes 
with the density-level. 

The morphology of the superclusters is 
quantified in terms of a s et of Shapefinders 
iSahni. Sathvaprakash fc ShandarinI 119981: ISheth et alJ 
The dimensional Shapefinders have dimensions of 
length and are given by 
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(9) 



For objects with fairly simple topology, the above Shapefind- 
ers give us a feel for the typical thickness (T), breadth (B) 
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Figure 7. Global Minkowski functionals are studied here w.r.t. FFm, which is defined as the overdense mass fraction. The MF- values 
are normalised to a volume of [100 h~^Mpc]^. In addition to the fact that the models can be distinguished, we also note that the geometry 
and topology of the excursion sets, when coupled with the overdense mass fraction, is markedly different in both the models at and 
beyond the onset of percolation. 



and length (L) of the object being studied. Since we have 
demonstrated that most of the superclusters are simply con- 
nected objects near or before the onset of percolation, this 
interpretation is valid for structures defined in this regime. 
Further, a combination of (T,B,L) can be used to define 
Planarity (P) and filamentarity (F) of the object. 



B-T 
B + T' 



F = 



L-B 
L + B' 



(10) 



Individual superclusters in a given realization of, say, 
tCDM may follow a distribution of shapes and sizes, not 
necessarily coinciding with a similar list due to another real- 
ization. Hence, to make a statistical study of a given model, 
we average over the shapes of individual superclusters de- 



fined at a given threshold of density, and test whether the 
average morphology as a function of the density-level is well 
constrained. For this purpose, we use the volume and mass 
averaged Shapefinders defined as follows. 



Pv,M (Sth) 
Fv,M (Sth) 



(11) 



where the quantity Qi may refer to the mass or volume of 
the i"' supercluster defined at the level 5th, and the summa- 
tion is over all the objects. It is to be noted here that the 
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Figure 8. The largest cluster spans most of the sample- volume for FFv>0.3. Here we study the dependence of the mean curvature and 
genus of the largest cluster of rCDM and ACDM w.r.t. FFv (left panels) and FFm (right panels). Note the well-constrained behaviour of 
both the measures for rCDM, and the distinctly different trend observed for them in ACDM. Please refer to the text for the discussion. 



average planarity and filamentarity are mostly contributed 
by the structures with large volume and/or mass. There are 
a finite number of such objects before the onset of perco- 
lation. Hence, in this regime, we probe the average shape 
of these objects. However, once the percolation sets in, the 
average morphology mainly reflects the morphology of the 
largest object in the system. 

We evaluate Pv,m and Fv,m for all the 11 samples at 
the density levels corresponding to (5€[0,(5max]. The measure- 
ments due to 10 rCDM realizations are used for averaging 
and for evaluating the errors. Figure|U]shows Pv,ai and Fv.m 
studied as functions of FFy (left panels) and FFm (right 
panels) respectively. The solid lines refer to rCDM and the 



dashed lines refer to ACDM. The vertical bars stand for the 
1(7 errors. We observe that 

• Pv,M and Fv,M are well constrained statistics and 
hence, can be used to reliably assess the morphology of the 
LSS. 

• The average planarity of the LSS steadily increases as 
we lower the threshold of density. 

• The average fllamentarity initially increases, reaches a 
maximum where the percolation sets in, but drops to sub- 
zero values after FFv~0.18. 

This implies that the square of the area of the excursion 
set is always greater than thrice the product of the volume 
and curvature. Further, is superseded by 47r(G-f 1)S at 
FFv>0.18. A steady drop in the curvature is responsible 
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Figure 9. The dimcnsionlcss Shapcfindcrs (P,F) arc evaluated for all the clusters at a given density level <5th- These Shapefinders are 
volume-averaged and are studied w.r.t. FFv (left panels). They are also inass-averaged and are studied w.r.t. FFm (right panels). The 
average shape of the cosinie web assumes increasingly planar morphology as is lowered. The corresponding negative filamentarity is 
due to the presence of large voids which lead to the development of apijreciable negative curvature in the system. The average filamentarity 
in the system is highest near the onset of percolation, thus making it ideal for the study of structures which are evident in the visual 
impression of the cosmic web. The average morphology revealed in terms of Pv,M and Fv,m is extremely well constrained for a given 
model and could be a useful diagnostic for discriminating the models of structure formation. 



for this behaviour. We note that the system also tends to 
become increasingly spongy. Taken together, these observa- 
tions lead to the following picture: as the density level is 
brought down, the overdense structures connect up across 
the sample-volume. This at a time loads to large scale fil- 
aments and also to the birth of large voids in the system, 
signifying the onset of percolation. As revealed in the top 
panels, the filamentarity of the system is highest at the 
time. With a further drop in 5, the filaments surrounding 
the "voids" grow thicker and locally these assume a slab- 
like structure. Such an excursion set, at a time has large 



planarity, but also fias largo negative curvature due to being 
surrounded by large, almost spherical voids. This negative 
curvature is responsible for the large negative filamentarity 

19 

We notice that both the statistics are extremely well- 



constrained. Further, 



This 



It was anticipated that the Shapefinders (P,F)€[0,1]. Our 
present exercise shows that excursion sets corresponding to <5<1 
show negative filamentarity. These objects are necessarily domi- 
nated by large negative curvature. 
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feature can be used to distinguish these two models, and is 
an instance of the fact that morphology of LSS is an equally 
robust statistic as MFs to compare LSS due to rival models 
of structure formation. 

We finally note that the superclusters defined according 
to the conventional definition {Ssc €[1,2]) do correspond to 
the percolating structures in the system, and indeed exhibit 
the highest filamentarity (the top panels of Figure (SJ. This 
proves that the Shapefinder statistics conform to our visual 
impression while quantifying the LSS. In the next subsec- 
tion, we shall try to quantitatively asses the sizes of these 
superclusters and shall demonstrate that these measures are 
sensitive to the models being investigated, enabling us to 
distinguish the two models successfully. 

6.4 Cumulative Probability Functions 

In this subsection we make a quantitative investigation of 
the sizes of the superclusters and study the variation of these 
measures as functions of the density contrast. Our main fo- 
cus is to make a statistical study of these measures and to 
develop a methodology to handle various statistics due to 
many realizations of a given model. We also compare the 
sizes of superclusters in rCDM with those in ACDM and 
find that these are sensitive to the model being investigated 
and can be used to distinguish the two models. 

For reasons explained earlier, we find it advisable to 
carry out the exercise over a set of levels of density. We 
study the sizes of superclusters starting from the highest 
density level ((5~3.5) down to the percolation threshold. 

Since the sizes of individual structures may vary from 
one realization to another, it is not meaningful to prepare 
a list of, say, top N superclusters in a given realization 
and check its statistical significance over other realizations. 
There may not be any unique correspondence of the rank 
of the supercluster in the list and its size. This is mainly 
because at the density-levels above or near the percolation 
threshold, there may be many structures with comparable 
sizes and masses. So the sorted array of structures in this 
range of density-levels will show a large amount of Pois- 
son fluctuations in the size of, say, the i"' largest cluster. 
Hence, a more meaningful approach is to study the cumula- 
tive probability functions of the Shapefinders so as to know 
how many superclusters in a given ensemble of them, have 
their sizes greater than a given triplet (T,B,L). 

As noted in Section 4, all the 10 samples of rCDM 
being studied here have comparable number of clusters at 
any given threshold of density, and that these are only 
marginally smaller in number than in ACDM. Further, we 
found that the global geometry and topology of the rCDM 
density fields is extremely well constrained. Hence, it is jus- 
tified to construct cumulative probability functions (here- 
after, CPFs) of various quantities using an ensemble which 
contains the clusters due to all the 10 realizations of rCDM 
The CPF of a given quantity Q is a fraction of the to- 
tal number of structures having the value of Q greater than 
a chosen value q. In our case, the total number of clusters 

However, prior to this, we have checked that the CPFs due to 
individual realizations of rCDM show trends similar to CPFs due 
to the cumulative ensemble. 



amounts to the sum over all the realizations of rCDM and 
is close to a few thousand. It turns out that it is beneficial 
to work with such large number of clusters, so that we have 
sufficient number of structures which are otherwise rare in 
a given sample. By finding these in sufficient number in a 
bigger ensemble, we can put an upper limit on the sizes 
of superclusters at any given threshold of density, without 
hampering our conclusions by Poisson noise. 

To begin with, we study the CPFs of the dimensional 
Shapefinders, the so called. Thickness (T), Breadth (B) and 
Length (L) of the superclusters with FFy as a parameter. 
Figure IIUI shows the results. We notice that the CPFs fol- 
low a specific pattern of evolution. While a typical range of 
thickness and breadth of the structures is fixed (T £(0,13) 
h"^Mpc and B G (0,17) h"^Mpc), the fraction of clusters 
with a given thickness and breadth increases with increasing 
FFy. Thus, the structures become thicker and broader as we 
lower the threshold, as anticipated. The CPF of length, on 
the other hand, extends to high L— values upon increasing 
FFv. This implies that the longer structures become statis- 
tically more significant with the decrease of 5, while their 
thickness and breadth increases only marginally. Thus, as 
5— >5pcrc, the filamentarity of the individual superclusters 
should increase. At the highest thresholds, however, the su- 
perclusters might look only marginally filamentary (F~0.25, 
P~0.1) (see. Figure]^, and in realistic situation, may give 
us an impression of being more like ribbons. . We infer that 
the top 10 superclusters in rCDM should have their length 
^40 h~^Mpc at 5^5max, and these should grow to become 
as long as L^90h~^Mpc at the onset of percolation. The 
longest superclusters at percolation, which are necessarily 
rare (say, a few in 1000), might be as long as 150h~^Mpc. 

To summarise, the supercluster-sizes are dependent on 
the threshold of density at which the superclusters are de- 
fined. Based on the study of CPFs as functions of 5th we infer 
that the thickness and breadth of the structures are compa- 
rable to one another at all thresholds of density 5^(5pcrc, 
with breadth B always slightly higher than thickness T. 
We have noted that at percolation, the largest superclusters 
could be as long as 150 h^^Mpc, and the superclusters with 
length Lg [90,150] h~^Mpc are relatively more frequently to 
be found (say, with a few per cent probability in an ensemble 
of 1000 structures). 

We are limited by the availability of only one realiza- 
tion of ACDM. Hence similar conclusions cannot be reliably 
drawn for this model, but a study of CPFs of ACDM super- 
clusters does reveal similar trends. 

In order to see whether the supercluster sizes and their 
morphology can be utilised to distinguish the two models, we 
studied the CPFs of T, B and L due to ACDM and rCDM 
together. We find that near to the onset of percolation, the 
two dimensional Shapefinders thickness and breadth, and 
the planarity of the superclusters do give a detectable signal. 
Figure 1111 shows the results. The thick solid lines refer to 

A more detailed study of the morphology in this range of 
5 needs to be carried out to understand whether the super- 
structures with significant planarity are statistically significant 
or not. This question has an observational bearing because some 
of the super-structures in our local neighbourhood appear planar. 
For example, the Great Wall, is characterised by (T,B,L)~(10, 50, 
120) h-^Mpc. See, Fairall 1998 and Martinez & Saar 2002. 
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Figure 10. The CPFs of dimensional Shapefinders are studied as functions of the density level for tCDM. The bottom panel shows 
the associated volume filling fractions. The percolation curves are shown as bold, solid lines. We notice that the clusters grow thicker, 
broader and longer as we increase the volume fraction (at any given triplet (T,B,L), there are more number of them thicker than T, 
broader than B and longer than L as we decrease 5). 



rCDM. The la error-bars are computed assuming a Poisson 
statistics and go as the square- root of the number of clusters 
available at a given value of (T,B,P). It has been checked 
that the CPFs due to individual realizations fall within la 
of the CPFs due to the total ensemble in most cases. Thus 
our decision to use the total ensemble (cumulative of all the 
clusters and due to all the realizations) and of assuming 
Poisson errors is justified. To make the comparison ideal, 
we choose to incorporate Poisson errors on the CPFs of the 
only available realization of ACDM as well. As shown in 
the figure, the rCDM CPFs of thickness and breadth are 
systematically higher than those of ACDM at FFv=9.1, 10.4 
and 11.7 per cent. Thus, we are able to establish that the 
rCDM superclusters with a given thickness are statistically 
larger in number than in ACDM. Assuming similar number 
of clusters in the two models, we could infer that the rCDM 
superclusters are thicker than those in ACDM. This effect 
is pronounced enough to be detected with a set of limited 
number of realizations. The topmost panels further reveal 



that CPF^(P)<CPF^(P). We conclude from here that the 
thickness, breadth and planarity of superclusters follow well- 
defined distributions which are special to the model being 
studied. 

In order to test whether we can constrain the length 
of superclusters of a given model, we studied the CPF for 
length L, following identical procedure. Figure [T?l shows the 
results, where the CPF(L) have been evaluated at FFv=7.8 
per cent (at the onset of percolation) and at FFv=11.7 per 
cent (after the percolation). We notice that at the onset 
of percolation, the length of the longest rCDM superclus- 
ters could be as large as 90 h~'^Mpc, whereas their ACDM 
counterpart structures, which exhibit same degree of sta- 
tistical significance, are relatively shorter with Lmax ^ 55 
h^^Mpc. We conclude that the rCDM superclusters tend to 
be statistically longer than their ACDM counterpart struc- 
tures. Thus, we can also utilise the length of the superclus- 
ters to compare and distinguish the two models. After the 
percolation, the longest structure dominates the survey vol- 
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Figure 11. The cumulative probability functions for thickness (T), breadth (B) and planarity (P) are shown for both the models at 
three thresholds of density. The fraction of clusters with a given thickness or breadth in rCDM are always greater than or equal to the 
fraction of clusters in ACDM. Further CPF^(P)>CPF^ (P). These Shapefinders can hence, be used to distinguish these models. 



ume and enters only as a Poisson fluctuation in the CPF 
of length. The remaining structures are evidently shorter in 
length compared to their longer progenitors at percolation, 
so that CPF(L) are confined to lower values of L. This effect 
is evident in the right panel of Figure lT^ In both the panels, 
shown as dashed line is the CPF(L) due to the first realiza- 
tion of tCDM which shares its set of initial random numbers 
with the ACDM simulation. We note the sharp tendency of 
tCDM structures to be longer than those of ACDM, an effect 
which successfully enables us to distinguish the two models 
of structure formation. 



7 DISCUSSION AND CONCLUSIONS 

This paper provides a theoretical framework and a method- 
ology to analyse the galaxy-distribution revealed by large 
3- dimensional redshift surveys. The exercise is in the wake 
of preparing to tackle the redshift surveys like 2dFGRS and 
SDSS, and consists in studying the geometry, topology and 



morphology of the LSS as revealed by the volume limited 
samples derived from 10 realizations of mock tCDM— based 
SDSS catalogue. These volume limited samples are prepared 
so as to refer to the same physical volume, and exhibit simi- 
lar number density of galaxies. The density fields used in our 
calculations are prepared after smoothing all the samples at 
6 h^^Mpc. 

We use SURFGEN (Sheth et al.2003) for calculating 
MPs, and find MPs as well as the derived morphological 
statistics, the Shapefinders, to be well-constrained statistics, 
useful to probe the contribution of the higher order correla- 
tion functions on the clustering of the galaxy-distribution. 
For this purpose, we compare the measured MPs due to 
10 realizations of rCDM with those due to a realization 
of ACDM. All the 11 realizations, by construction, repro- 
duce the observed two— point correlation function measured 
using the APM catalogue, and show similar clustering am- 
plitude. Some of the important results of our analysis are 
summarised below: 
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Figure 12. Shown here are the cumulative probabihty functions of length L for clusters of rCDM and ACDM at two thresholds of 
density corresponding to FFv=7.8 per cent (just before the onset of percolation) and FFv=ll-4 per cent (after the onset of percolation). 
We find the CPFs of the two models to be distinctly different at longer length-scales. The dashed line refers to the CPFs due to the first 
realization of tCDM which shares the initial set of random numbers with ACDM. The largest superclusters of tCDM are systematically 
larger than those due to ACDM. 



• We show that the global MFs (defined as the sum of par- 
tial MPs over the set of all the clusters at a given threshold 
of density) of rCDM show systematically lower amplitude 
compared to those due to ACDM, an elTect which enables us 
to distinguish the two models, and can be attributed to the 
nonzero higher order correlation functions. Considering the 
fact that higher order correlation functions are cumbersome 
to estimate and offer little insight in their interpretation, this 
can be termed as a success of MFs as efficient quantifying 
statistics of LSS. 

• The Virgo simulations of dark matter showed higher 
amplitude of MFs for rCDM compared to ACDM (see Sheth 
et al.2003). The cosmological parameters used by Virgo 
group match with those employed by CHWF except for the 
shape parameter F, which is 0.21 in the former case and 0.25 
in the latter. We assume that this minute difference in F will 
not alter the above trend so long as both the models start 
with the same power spectrum. The galaxy distributions due 
to both the models were analysed in this paper, and we found 
a reversal in the above trend: The MF-amplitudes for rCDM 
were found to be smaller than ACDM. 

We note that the relative smallness of the amplitudes 
of MFs reflects the higher degree of phase correlations in 
the matter-distribution. Taken together with this fact, our 
above observation points out that scale-dependent biasing 
could lead to an apparent phase-mismatch between the 
galajcy-distribution and the underlying matter-distribution. 
Thus, there may be larger degree of phase-correlations in 
the matter- distribution of ACDM compared to rCDM, but 
the biased galaxy-distribution of ACDM would appear to 



exhibit smaller phase-correlations compared to rCDM. We 
note however, that these results are afi'ected by the con- 
straint that the bias is fixed so as to reproduce the observed 
2— point correlation function (see Sectio n 2). 

An analysis of 2dFGRS performed bv lLahav et al.l ll2002l) 
concluded that the average bias on scales 0.02 fc 0.15 h 
Mpc~^ is indeed ~1 In light of the above fact, it will be 
interesting to see how the geometry, topology and morphol- 
ogy of the LSS change in the presence of scale-independent, 
constant bias. From our analysis, however, it is clear that 
bias in various structure formation scenarios can play im- 
portant role, and comparing the dark matter distributions 
with the SDSS and/or 2dFGRS galaxy distributions is not 
an advisable way of testing the model(s) of structure forma- 
tion against the observations. A realistic treatment of bias 
en route the proper treatment of the physics of galaxy for- 
mation in preparing a mock catalogue, is desired in any such 
exercise. 

• Since we had access to 10 realizations of rCDM, we 
made a statistical study of the shapes and sizes of the su- 
perclusters occurring in this model. We found the thickness 
and breadth of the superclusters to increase only marginally 
over a range of S (T, BG[l,17]h"iMpc; TsCB). However, the 
length of the largest superclusters increases monotonically 
from ~40 h"^Mpc to ~90 h"^Mpc as S^Spcic- Near the per- 
colation, the rCDM superclusters with length Js90 h~^Mpc 
are ~ 1 per cent probable to be found (say, the top 10 in 

However, it should be noted that the cosmological model was 
assumed to be ACDM. 
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an ensemble of 1000 structures). However, the longest (and 
necessarily the most massive) superclusters could be as long 
as 150 h~^Mpc and are rare to be found, with less than 1 
per cent probability in a volume as big as that covered in 
this analysis. These findin gs are to be confron ted with re- 
cent observational work bv lBrand et al.l JioO^) who report 
rare structures of the size of ~100 Mpc being traced by ra- 
dio galaxies. The dominant morphology of the superclusters 
is prolate, or ribbon-like at the highest threshold of density, 
and it evolves to become more filamentary as the threshold 
is lowered toward the percolation threshold. 

• We further found that the CPFs of thickness and 
breadth near the percolation threshold show distinctly 
higher amplitude in rCDM compared to ACDM. Even pla- 
narity of the two models shows appreciably different distri- 
bution for the two models. These Shapefinders, thus, could 
be good discriminants of the structure formation scenarios. 

• A study similar to the above, repeated for the length 
of the superclusters reveals that the longest superclusters of 
rCDM are likely to be longer than their ACDM counterpart 
structures at the same level of significance. In particular, 
whereas the longest superclusters of rCDM could be as long 
as 90 h"^Mpc, the ACDM superclusters have L>55 h"^Mpc. 

To summarise, we have established that the Minkowski 
Functionals estimated using SURFGEN are well constrained 
statistics and can be reliably used to quantify the geometry 
and topology of the LSS in various scenarios of structure 
formation. We have used these to distinguish two rival mod- 
els of structure formation, namely ACDM and rCDM, and 
have explored a set of related measures and morphologi- 
cal statistics like Shapefinders to make such a comparison 
more robust. We have further measured and reported the 
sizes and shapes of the superclusters in rCDM, with a scope 
of extending similar methodology to other models with the 
availability of larger number of realizations for them. 

Our present exercise does have its limitations, which are 
mentioned below. 

• The efi'ect of redshift space distortions has not been 
incorporated in the present exercise. We have throughout 
worked in the redshift space and have not quantified the 
fingers-of-God effect. Such a treatment is beyond the scope 
of this paper and shall be taken up in a future work. 

• The mock catalogues studied in this paper are pre- 
pared by implementing ad-hoc biasing schemes. Although 
the methodology of analysis is clearly laid out in the present 
work, a more realistic comparison of the mock data with the 
observed LSS of the Universe should be carried out using 
the mock catalogues prepared by incorporating the physics 
of galaxy-formation in as much detail as possible. The semi- 
analytic methods, the SPH simulations of galaxy formation 
and the Conditional Luminosity Function based approach 
should be able to deliver such mock catalogues in near fu- 
ture. 

We hope to address some of the issues mentioned above and 
extend this work to a comparison of the mock catalogues 
with the LSS revealed by 2dFGRS and SDSS in near future. 
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